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ABSTRACT 



We consider some models of f{R) gravity that can be used to describe, in a 
suitable weak-field limit, the gravitational field of the Sun. Using a perturba- 
tive approach, we focus on the impact that the modifications of the gravitational 
field, due to the non-linearity of the gravity Lagrangian, have on the Solar System 
dynamics. We compare the theoretical predictions for the precession of the lon- 
gitude of the pericentre w oi a, test particle with the corrections to the standard 
Newtonian-Einsteinian precessions of the longitudes of perihelia of some planets 
of the Solar System recently estimated by E.V. Pitjeva by fitting large data sets 
with various versions of the EPM ephemerides. 

Subject headings: Experimental tests of gravitational theories; Modified theories of 
gravity; Celestial mechanics; Orbit determination and improvement; Ephemerides, 
almanacs, and calendars 
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Introduction 



General Relativity (GR) has passed with excellent results many observational tests: 
a satisfactory agreement comes both from S olar System tests and from binary pulsars 



observations. As a matter of fact (see e.g. (jWill 



20061 )). the current values of the PPN 



parameters are in agreement with GR predictions and, consequently, Einstein's theory is 
the classical theory of gravitational interactions accepted nowadays. 

However, observations seem to question the general relativistic model of gravitational 
interactions on large scales. On the one ha nd, the data coming from the rotation curves of 



spiral galaxies flBinney and Tremaine 



19871 ) cannot be explained on the basis of Newtonian 
gravity or GR: the existence of a peculiar form of matter is postulated to reconcile the 
theoretical model with observations, i.e. dark matter, which is supposed to be a cold and 
pressureless medium, whose distribution is that of a spherical halo around the galaxies. 
Furthermore, dark matter can explain the mass discrepancy in galactic clusters (jClowe et al. 



20061 ). On the other hand, a lot of observations, such as the light curv es of the type la 



supernovae and the cosmic microwave bac kground (CMB) experiments f lRiess et al. 



Perlmutter et al. 



1998 



1999 



Bennet et al. 



20031 ). firmly state that our Universe is now undergoing 
a phase of accelerated expansion. Actually, the present acceleration of the Universe cannot 
be explained, within GR, unless the existence of a cosmic fluid having exotic properties is 
postulate, i.e. dark energy or introducing a cosmologic al constant which, in tur n, brings 



about other problems, concerning its nature and origin (jPeebles and Ratra 



20031 ). 



The main problem one has to face with dark matter and dark energy (or the 
cosmological constant) is understanding their nature, since they are introduced as ad hoc 
gravity sources in GR or its weak-field limit, Newtonian gravity. 

In order to explain the observations another possibility exists: the query for dark 
matter and dark energy points out the failure of GR (and its approximation, Newtonian 
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gravity) to deal with gravitational interaction at galactic, intergalactic and cosmological 
scales. The latter viewpoint led to the introduction of various modified gravity models. 

In this paper, we are concerned with the so called f{R) theories of gravity, where 



the gravitational Lagrangian 



(jCapozziello and Francaviglia 



depends on a function f o. 



2007 



Sotiriou and Faraoni 



the scalar curvature R (see 



20081 ) and references therein). 



These theories are also referred to as "extended theories of gravity", since they naturally 
generalize GR: in fact, when f{R) = R the action reduces to the usual Einstein- Hilbert 
action, and Einstein's theory is obtained. These theories can be studied in the metric 
formalism, where the action is varied with respect to metric tensor, and in the Palatini 
formalism, where the action is varied with respect to the metric and the affine connection, 
which are supposed to be independent from one another (actually, there is also the 
possibility that the matter part of the action depends on the affine connection, and is 
then varied with respect to it: this is the so-called metric- affine formalism, but we are 
not concerned with this approach in this paper). In general, the two approaches are not 
equivalent: the solutio ns of the Palati ni field equations are a subset on solutions of the 



metric field equations (IMagnano 



199J). 



Actually, f{R) theories provide cosmo logically viable models, where both the 



infiat i on pha se and the accelerated expansion are reproduced (see (INojiri and Odintsov 



2007 



2008a 



b|) and references therein). Furthermore, they ha ve been used to explain 



the rotation curves of galaxies wi thout need for dark matter (jCapozziello et al 



Frigerio Martins and Salucci 



2003) 



2007al : 



However, because of the excellent agreement of GR with Solar System and binary 
pulsar observations, every theory that aims at explaining galaxies dynamics and the 
accelerated expansion of the Universe, should reproduce GR at the Solar System scale, 
i.e. in a suitable weak-field limit. In other words, also for f{R) theories the constraint 
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holds to have correct Newtonian and post-Newtonian hmits. This issue has been hvely 
debated in the recent hterature, where different approaches to the problem have been taken 



into account, both in the Palatini and metric formal 
found in the recent review by 



ism. A thorough discussion can be 



Sotiriou and Faraonil (120081 ). In summary, with res pect to 



the w eak-field tests and, more in general, the non cosmological solutions (see e.g. (IFaraoni 
20081 )). it seems that there are difficulties in considering Palatini f{R) gravity as a viable 
theory because the Cauchy problem is ill-posed and, furthermore, curvature singularities 
arise when dialing with simple stellar models; as for metric f{R) gravity, there are models 
that are in agreement with the weak-field tests, but it seems that curvature singularities 
exist, in this case, for compact relativistic stars. 

Without going into the details of this interesting debate, in this paper we want to 
test some models of f{R) gravity that can be used to describe the gravitational field of 
the Sun. In particular, we are going to examine the impact that the modifications of the 
gravitational field of GR have on the Solar System dynamics. We apply a perturbative 
approach to compare the /(i?) -induced secular precession of the longitude of the pericentre 
Ct7 of a test particle with the latest determinations of the corrections to the usual perihelion 



precessions c oming from fits o 



ephemerides (jPitjeva 



2005a 



• huge 



2006 



planet ary data sets with various versions of the EPM 



2008a 



a) 



The paper is organized as follows: in Section [2] we briefly review the theoretical 
formalism of f{R) gravity, both in the metric and Palatini approach, then, in Section [3] we 
outline a general approach to the perturbations of the gravitational field of GR, due to the 
non-linearity of the gravity Lagrangian. In Section |4] we compare the theoretical predictions 
with the observations. Finally, discussion and conclusions are in Section [51 
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2. The field equations of f{R) gravity 
In this Section, we introduce the field equations of f(R) gravity. We shall consider 



both the metric and the Palatin i approach (see, e.g., (jCapozziello and Francaviglia 



and (jSotiriou and Faraoni 



20081)) 



20071 ) 



The equations of motion of f{R) extended theories of gravity can be obtained by a 
variational principle, starting from the action: 

A = Agrav + ^mat = j iVd f (R) + 2x^mat(^, V^/^)] d^'x. (1) 

The gravitational part of the Lagrangian is represented by a function f{R) of the scalar 
curvature R. The total Lagrangian contains also a first order matter part L^g^t, functionally 
depending on matter fields \E', together with their first derivatives, equipped with a 
gravitational coupling constant x = the metric formalism, F is supposed to be the 

Levi-Civita connection of g and, consequently, the scalar curvature R has to be intended 
as i? = R{g) = g^^Raisig)- On the contrary, in the Palatini formalism the metric g and 
the afiine connection F are supposed to be independent, so that the scalar curvature R has 
to be intended as i? = R{g,T) = g^f^RapiV), where R^yiV) is the Ricci-like tensor of the 
connection F. 

In the metric formalism the action ([T]) is varied with respect to the metric g, and one 
obtains the following field equations 

f\R)R,. - \f{R)g,u - iS.Vu - g,u^) f'{R) = ^ T^., (2) 

where f{R) = df{R)/dR, and T^" = —-^^j^ is the standard minimally coupled matter 
energy-momentum tensor. The contraction of the field equations ([2]) with the metric tensor 
leads to the scalar equation 

3n/'(i?) + nR)R - 2f{R) = ^T, (3) 



- 7- 



where T is the trace of the energy-momentum tensor. Eq. ([3]) is a differential equation for 
the scalar curvature R. 

In the Palatini formalism, by independent variations with respect to the metric g and 
the connection F, we obtain the following equations of motion: 

/'(i?)%.)(r)-i/(i?)^7,. = ^T,., (4) 

Vl[^gnR)gn = 0, (5) 
where V'" means covariant derivative with respect to the connection F. Actually, it is 



possible to show ([Ferraris et al 



1993 



19941 ) that the manifold M, which is the model of 
the space-time, can be a posteriori endowed with a bi-metric structure (M , g, h) equivalent 
to the original metric-affine structure (M, (7, F), where F is assumed to be the Levi-Civita 
connection of h. The two metrics are conformally related by 

h^u = f'iR) gt^u- (6) 

The equation of motion (j4]) can be supplemented by the scalar-valued equation obtained by 
taking the contraction of (j4]) with the metric tensor: 

f\R)R-2fiR) = ^T. (7) 

Equation ([7]) is an algebraic equation for the scalar curvature R. 

In order to compare the predictions of f{R) gravity with Solar System data, we have 
to consider the solutions of the field equations ([2]) , (j4]) , ([5]) - supplemented by the constraints 
([3]), (IT!) - in vacuum, since tests are based on the observations of the dynamics of the planets 
in the gravitational field of the Sun. 

In particular, the vacuum field equations in the metric approach read 

f'{R)R,, - \f{R)g,. - (V,V, - g,,U) f\R) = 0, (8) 
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supplemented with the scalar equation 

3Df'{R) + f{R)R-2f{R)=0. (9) 
In the Palatini approach, the field equations become 

f'(R)R^^^^(T)-^f{R)g,, = 0, (10) 

'^l[^/^f'{R)9n = 0, (11) 
and they are supplemented by the scalar equation 

f'{R)R-2f{R)=0 (12) 

We want to point out some general features of the scalar equations iQ and f|T2l) . which can 
help to understand the differences between the vacuum solutions in the two formalisms. 



In Palatini f{R) gravity, the trace equation 



admits constant solutions R = Ci ( iFerraris et al. 



([T^ is an algebraic equation for R, which 



19931 ). and it is identically satisfied if f{R) 



is proportional to R^. As a consequence, it is easy to verify that (if f'{R) 7^ 0) the field 
equations become 

Rfiu = -^Rdfiu (13) 

which are the same as GR field equations with a cosmological constant. In other words, in 
the Palatini formalism, in vacuum, we can have only solutions that describe space-times 
with constant scalar curvature R. Summarizing, eq. f|T3l) suggests that all GR solutions 
with cosmological constant are solutions of vacuum Palatini field equations: the function 
f{R) determines the solutions of algebraic equation (fT2l) . 

In metric f{R) gravity the trace equation ([9]) is a differential equation for R: this 
means that, in general, it admits more solutions than the corresponding Palatini equation. 
In particular, we notice that if i? = constant we obtain the Palatini case: so for a given 
f{R) function, in vacuum, the solutions of the field equations of Palatini f{R) gravity are 
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a subset of the solutions of the field equations of metric f{R) gravity (IMagnano 



19941 ) 



howe ver, in metric f(R) gravity , vacuum solutions with variable R are allowed too (see, 



e.g., (IMultamaki and Vilja 



20061)) 



3. Corrections to the gravitational potential 

We have seen in the previous Section that, when f{R) ^ R, the field equations of f{R) 
gravity are different from those of GR. Thus, it is evident that the solutions of such modified 
field equations describing the gravitational field of a point-like mass (e.g. the Sun) contain 
corrections to the GR solutions, both at Newtonian and post-Newtonian level. However, 
these corrections have to be small enough not to contradict the known tests of GR. Thus, 
it is possible to treat them perturbatively to evaluate their impact on the dynamics of the 
Solar System planets. 

In this Section we want to outline the general procedure that we are going to apply to 
some solutions of j[R) gravity that can be used to describe the gravitational field of the 
Sun, in order to compare the predictions of these gravity models with the existing data. 

In general, we are going to deal with spherically symmetrical metrics, describing 
the space-time around a point-like mass M, which can be endowed with proper angular 
momentum J . The weak-field and slow-motion approximations of these metrics will be 
sufficient for our purposes. Generally speaking, this means that the deviations from GR 
will be linear in some parameters deriving from the specific j[K) gravity model. 

On using spherical isotropic coordinates, these metric have the general formic 

ds^ = A{r)df + B{r) {dr^ + r^d^^ + sin^ M^^) + 2C{r)sin^edtdip, (14) 



^If not otherwise stated, here and henceforth we use units such that G = c = 1. 
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where the angular momentum J is assumed to be perpendicular to the 6* = 7r/2 plane. 

The gravitational (scalar) potential $(r) is read from the A{r) function 

A(r) = l + 2$(r). (15) 

According to what stated before, we expect a gravitational potential in the form 

$(r) = $^(r) + A$(r), (16) 

where $^(r) = — ^ is the Newtonian potential of a point-like mass M, and A$(r) ^ $Af(r) 
is a correction vanishing for f{R) R. 



The C(r) function accounts for the presen ce of the so-called gravito-magnetic effects 



flRuggiero and Tartaglia 



2002 



Mashhoon 



20071 ) induced by the rotation of the source of the 



gravitational field. In GR C (r) is given by the suitable component of the gr ayito- magne tic 



vecto r potential of a gravito-magnetic dipole, i.e. A^^(? 



2J 



see e.g. (IMashhoon 



20071 )). As a consequence, we expect that the C{r) function has the form 



C(r) = Aj» + AA^(r) 



(17) 



where, again, ^A^p{r) -C A^p{r) is a correction vanishing for f{R) — >■ R. 



We can use the gravito-electromagnetic (IRuggiero and Tartaglia 



2002 



Mashhoon 



20071) 



formalism to describe the total perturbing acceleration felt by a test particle in the metric 

m 



W = -S^ -2vxB 



G 



where 



and 



rfA$(r) 
dr 



= V X A, A - 



r sin 9 



(19) 



(20) 
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Hence, given the perturbing acceleration f|T8|) . we can calcu 



ate it s effects on planetary 



Rovl (120051)). We may use 



motions within standard perturbative schemes (see, e.g. 
the Gauss equations for the variations of the elements, which enable us to study the 
perturbations of the Keplerian orbital elements due to a generic perturbing acceleration, 
whatever its physical origin is. The Gauss equations for the variations of the semi-major 
axis a, the eccentricity e, the inclination i, the longitude of the ascending node Q, the 
argument of pericentre u and th e mean anomaly of a test particle in the gravitational 



field of a body M are (IRoy 



20051 ) 



da 
Itt 
de 
dt 
di 
It 

dn 

~dt 

doj 

m 

dM 
"df 



n 



eWr sin v + Wr{- 
r 



na 



Wr sin V + Wr 



1 / r 

cos V -\ — 1 

e V a 



_ ^ ^ (-) cos(w + v), 

——^-== (-) sm{u + v), 



cos z 



.dQ v^r^ 



dt 



nae 



2 /r 
n- — Wr{- 

na \a 



/ duj 
\dl 





r\ 


— smw 




PJ 


dn\ 









(21) 
(22) 
(23) 
(24) 

(25) 
(26) 



in which n = 27t/P is the mean motioiu, P is the test particle's orbital period, v is 
the true anomaly counted from the pericentre, p = a(l — e^) is the semilatus rectum of the 
Keplerian ellipse, Wr, Wr, Wy are the radial, transverse (in-plane components) and normal 
(out-of-plane component) projections of the perturbing acceleration W , respectively, on the 
orthonormal frame {f, f, €>} comoving with the particle. 

For our purposes it is useful to consider the longitude of the pericenter zu = u + cos i Q. 



'For an unperturbed Keplerian ellipse it is n = ^GM/a\ 



- 12 - 



The Gauss equation for its variation under the action of an entirely radial perturbing 
acceleration Wr is 



dw 
It 



-Wr cosv. 



nae 



(27) 



After being evaluated onto the unperturbed Keplerian ellipse, the acceleration (fTSjl must be 
inserted into eq. fl27|) : then, the average over one orbital period P must be performed. To 
this end it is useful also to recall the following relations where also the eccentric anomaly E 
is used 

r = a(l — e cos-E), 



dt 



(1— e cos E) 



dE, 



(28) 



cos w 



smf 



cos E—e 
1— e cos E ' 



sinE\/l—e'^ 
1—e cos E 



In fact, what we aim at is evaluating the perturbations induced on the longitudes of 
the perihelia by the corrections to the gravitational field due to f{R) gravity, in order to 
compare them with the latest observational determinations. The astronomer E.V. Pitjeva 
(Institute of Applied Astronomy, Russian Academy of Sciences, St. Petersburg) processed 
almost one century of data of different ty pes for the major bodi es of the Solar System to 



improve the EPM planetary ephemerides fjPitjeva 



2005al . 



2008a 



b|). Among other things. 



she simultaneou sly estimated c orrections to the secular rat es of the longit udes of perihelia 



w of the inner (jPitjeva 



2008 



b|) planets of 



2005bl ) and of some of the outer (iPitjeval 
the Solar System as fit-for parameters of global solutions in which she contrasted, in a 
least-square way, the observations to their predicted values computed with a complete set of 
dynamical force models including all the known Newtonian (solar quadrupole mass moment 
J2, A^— body interactions with the major planets, 301 biggest asteroids, massive ring of the 
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small asteroids, 20 largest trans-Neptunian objects and massive ring for the other ones) 
and Einsteiniaijf] features of motion. As a consequence, any force that is not present in 
Newtonian gravity or GR is, in principle, accounted for by the estimated corrections to the 
usual apsidal precessions. For the sake of completeness, we reproduce in tables [1] and [2] the 
estimated perihelia extra-precessions for inner and outer planets, respectively. 

What we want to do is to see whether the estimated perihelia extra-precessions are 
compatible with the perturbations of the gravitational field deriving from the non-linearity 



Now, let us briefly outline how we are going to put f{R) gravity on the test. In 
general a correction to the gravitational field of GR due to the non linearity of the gravity 
Lagrangian, in the weak-field and slow motion approximation, can be parameterized in 
terms of a parameter k, where k ^ as far as f{R) R- In other words, k is a measure 
of the non-linearity of the Lagrangian. Let V{f{R)) be the prediction of a certain effect 
induced by these modified gravity models, e.g. the secular precession of the perihelion of a 
planet: for all the f{R) models that we are going to consider below, it turns out that 



where g is a. function of the system's orbital parameters a (semi-major axis) and e 
(eccentricity); such g is a peculiar consequence of the f{R) gravity model. Now, let us take 
the ratio of V{f{R)) for two different systems A and B, e.g. two Solar System's planets: 
VA^f (R)) fPsif (R)) = gx/QB- The model's parameter n has now been canceled, but we 
still have a prediction that retains a peculiar signature of that model, i.e. gx/gB- Of course, 
such a prediction is valid if we assume k is not zero, which is just the case both theoretically 
(only if f{R) = R then k = 0) and observationally because k is usually determined by other 



of /(i?). 




(29) 



■^The general relativistic gravito-magnetic Lense-Thirring force has not yet been modeled. 
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able 1: Inner planets. First row: estimated perihelion extra-precessions, from Table 3 of 



( jPitjeva 



2005bl ). The quoted errors are not the formal ones but are realistic. The units 



are arc-seconds per century (" cy ^). Secon d row: semi-rn aior axes, in Astronomical Units 



(AU). Their formal errors are in Table IV of (jPitjeva 



2005al ). in m. Third row: eccentricities. 



Fourth row: 


orbital periods in years. 








Mercury 


Earth 


Mars 


(^) (" cy- 


1) -0.0036 ±0.0050 - 


-0.0002 ±0.0004 


0.0001 ±0.0005 


a (AU) 


0.387 


1.000 


1.523 


e 


0.2056 


0.0167 


0.0934 


P (yr) 


0.24 


1.00 


1.88 



Table 2: Outer planets. First row: estimated perihelion extra-precessions (iPitjeval 120061 ). 
The quoted uncertainties are the formal, statistical errors re-scaled by a factor 10 in order 
to get the realistic ones. The units are arc-seconds per century (" cy~^). Second ro w: semi- 



ma i or axes, in Astronomical Units (AU). Their formal errors are in Table IV of (jPitjeva 



2005al ). in m. Third row: eccentricities. Fourth row: orbital periods in years. 



Jupiter 



Saturn 



Uranus 



(ro) (" cy"i) 0.0062 ±0.036 -0.92 ± 2.9 0.57 ±13.0 



a (AU) 
e 

P(yr) 



5.203 
0.0483 
11.86 



9.537 
0.0541 
29.45 



19.191 
0.0471 
84.07 
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independent long-range astrophysical/cosmological observations. Otherwise, one would 
have the meaningless prediction 0/0. The case k, = (or k < k, i.e. when k, is negligibly 
small) can be, instead, usually tested by taking one perihelion precession at a time. If 
we have observational determinations O for A and B of the effect considered above such 
that they are affected alsqj by the f{R) gravity model (it is just the case for the purely 
phenomenologically estimated corrections to the standard Newton-Einstein perihelion 
precessions, since any f{R) gravity model has not been included in the dynamical force 



models of the ephemerides adj usted to the p 



estimation process by Pitjeva (jPitjeva 



2005a 



anetary data in the least-square parameters' 
we can construct Oa/Ob and compare 



it with the prediction for it by /(-R), i.e. with qa/qb- Note that 50/0 > 1 only means 
that O is compatible with zero, being possible a nonzero value smaller than 60. Thus, it is 
perfectly meaningful to construct Oa/Ob- Its uncertainty will be conservatively evaluated 
as \1/Ob\SOa + \0a/0^\60b- As a result, Oa/Ob will be compatible with zero. Now, the 
question is: Is it the same for qa/qb as well? If yes, i.e. if 



Oa _ VA{f{R)) 
Ob VB{f{R)) 



(30) 



within the errors, or, equivalently, if 



Oa VA{f{R)) 







(31) 



Ob VB{f{R)) 

within the errors, the f{R) gravity model examined can still be considered compatible with 
the data, otherwise it is seriously challenged. 

In next Section some solutions of f{R) gravity that can be used to describe the 
gravitational field of the Sun will be tested according to the procedure that we have just 



'^If they are differential quantities constructed by contrasting observations to predictions 
obtained by analytical force models of canonical Newtonian/Einsteinian effects, O are, in 
principle, affected also by the mis-modeling in them. 
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described. 



4. f{R) weak-field solutions 

In this Section we introduce some solutions of f{R) gravity that have been used in the 
literature to the describe the weak gravitational field, and that can be considered as suitable 
models of the gravitational field of the Sun. We consider these modified gravitational fields 
and, within the perturbative scheme outlined above, compare the theoretical predictions 
with the estimated extra-precessions of the planetary perihelia. 



4.1. Power law corrections 



Capozziello et al. 



(l2007a|), in the 



Starting from a Lagrangian of the form f{R) = foR^' 
metric approach, look for solutions describing the gravitational field of a point-like source, 
in order to reproduce the galaxies rotation curves without need for dark matter. As a 
result, they obtain the following power-law form for the gravitational potential: 



M 
r 



1 



(32) 



The deviation from the Newtonian potential is parameterized by a power law, with two 
free parameters (3 and r^. In particular, (3 is related to n, i.e. the exponential of the Ricci 
scalar in f{R) = /oi?". 



In this case, the correction to the gravitational potential is 



M frV 

A$(r) = , 

r / 



(33) 
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which clearly leads to the radial acceleration 

Wr 



It yields the following perihelion precession (jlorio and Ruggierdl2008al ) 



with G(eA; /?)/G(eB; /9) ~ 1 for all the planets of the Solar System. 



Capozziello et al 



(34) 



(35) 



(j2007al ) find [3 = 0.817 from a successful fit of several galactic 
rotation curves with no dark matter; it is ruled out by comparing for seve ral pairs of planets 



Ari7A/Ari7B to Pa/^b obtained from eq. fl5^ fllorio and Ruggiero 



2008af ) 



4.2. Schwarzschild-de Sitter-like corrections 



The field equations (HHSl) and the stru ctural equation (j7|) , in t he Pa 



the s pherically symmetrical solution (see (lAUemandi et al. 



20071)): 



atini formalism have 



20051 ) and (IRuggiero and lorio 



1 



2M k 



1 + 



2M k 



{dr^ + r''de'' + 



r sm 



'■) , (36) 



r 3 / \ r 6 

k = Ci/i, where R = Ci is any of the solutions of the structural equation ([7]). In this case, 
we may write 

, , M kr^ 
m = + 37 

In this case, the perturbing potential is 

kr'^ 



A$(r) 



and the induced the perturbing acceleration is 



6 



(38) 



Wr 



--kr. 



(39) 
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As an y other Hooke-type extra-acc eleration, eq. (139|) induces a secular perihelion 



precession ( iKerr et al. 



2003 



lorio 



2008af ) 



(w) oc — = k\ —. 
^ ' n V M 



(40) 



By using eq. (HOl) to construct Va/Vb for different pai rs of Solar S ystem's planets and 



comparing them to ArijA/AtijB yield a negative answer (jlorio 



2008al ). 



4.3. Logarithmic corrections 



Soboutil (120071 ) aims at determining a f{R) able to explain the rotation curves of the 
galaxies obtained. In particular, working in the metric approach, solutions with R variable 
with the radial coordinate r are obtained. In this context, the gravitational potential reads: 



$(i?) = -—+2 + 2 



(41) 



The pararn eter a can be related to Modified Newtonian Dynamics (MOND, see e.g. 



(IMilgrom 



19831 )) characteristic acceleration Aq 



The Logarithmic correction to the Newtonian gravitational potential assumes the form 



A$(r) = -7Mln 



(42) 



and leads to a perturbing radial acceleration 



7M 



(43) 



In particular, in order to agree with the potential ( HTl) . we must set 7 = —a/2, ro = 2M. 



This kind of acceleration has been treated by 
approach outlined here getting negative answers. 



lorio and Ruggierd (l2008al ) with the 
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4.4. Yukawa-like corrections 



the Palatini ( 


Hamitv and Barracc 


1993; 


Barraco et al. 


1996) 


(Pechlaner and Sex] 


1966; 


Stelle 


1978 




Capozziello et al. 





2007bl )) Yukawa-like corrections are obtained. They lead to a gravitational potential in the 



form 



M 
r 



1 + a exp 



A 



(44) 



The parameter a is related to the strength of the correction, while A is related to the 
range of the modified potential. 



expi-- 



The Yukawa correction to the Newtonian potential 

A$(r) = -- 

yields an entirely radial extra-acceleration 

Ma r\ f r 



(45) 



(46) 



By only assuming A ^ ae, i.e. Yukawa-type long-range modifications of gravity, it is 
possible to obtain useful appro ximated exp ressions for the induced perihelion precession 



which, in turn, allow to obtain (lorio 



2007ah 



In 



I as Acta \ 
a\ AiiJB ) 



for the range and 



2A^Aro /a 
a = — , exp 



(47) 



(48) 



VMa VA 

for the strength. By using A = Earth, B = Mercury in eq. fl47p one gets A = 0.182 ± 0.18 3 
AU; su ch a value for A, the data of Venus and eq. (jUj) yield a = (—1 ± 4) x 10^^^ ( llorio 



2008bf ) 
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4.5. Gravito-magnetic effects 

We have shown that the vacuum solutions of General Relativity with a cosmological 
constant can be used in Palatini f{R) gravity. In particular, the Kerr-de Sitter solution, 



which describes a rotating black-hole in a space-time with a cosmo . 



(IDemianski 



1973 



Carter 



19731: 



Kerr et al. 



2003 



Kraniotis 



2004 



2005 



ogica l constant 



20071 ). can be used to 



investigate Gravito-magnetic effects in extended theories of gravity. 



In particular, the wea k-field and slow-motion approximation of the Kerr-de Sitter is 



(llorio and Ruggiero 



2008bh 



1 

+ 2 



de 



2M k 

r 3 

J f2M k 

— h -r 

M \ r 3 



2M k 



r 



6 



r^^ {dr^ + r^de^ + 



2 • 2/)>i2 

r sm ^''-^ 



) + 



2 ' -Mkr ] sin^ ed(j)dt. 



(49) 



We obtain the following expression for the perturbing gravito-magnetic potential 



AA^ = — i-r^ + -Mkr ] sin^O. 
^ M \3 6 



(50) 



Furthermore, the perturbing acceleration is 



W = -2vx B 



G 



(51) 



where the gravito-magnetic field B*^ is 



Q _ Jk - 5Jk 
^ "3M'^+ 12 



J + J • f f 



(52) 



The resulting orbital effects are (jlorio and Ruggiero 



2008bh 
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(d) 



0, 

0, 

0, 
Jk 



5M 



M 



3M V 2a 
2Jk cosi 



n + 



3M 
5Jk cosi 



3M 



5M 

4a 

M 
1 + — 
a 



(53) 
(54) 
(55) 
(56) 

(57) 

(58) 



In the calculation we have neglected terms of order 0{e^). 

By using the corrections Arij separately for each Solar System's planet one gets 

k < 10^29 



For all the Solar System's planets the perihelion rate can be satisfactorily approximated 



by 



2Jk cost 
3M ■ 



(59) 



Since cos^a/ coszb ~ 1 for every pair of planets A and B, in this case Va/V -r ~ 1; this 
possi bility is ruled out by ( 9a /(9r = A wa/ Awp,, as in the case of the DGP (IDvali et al 



2000l ) braneworld scenario (jlorio 



2007bf ) 



5. Discussion and Conclusions 

In this paper we have considered some solutions of f{R) gravity, both in the Palatini 
and metric formalism, that can be used to describe the weak gravitational field around 
the Sun. In particular, we have focused on the impact that the modifications of the GR 
gravitational field, due to the non-linearity of f{R), have on the Solar System dynamics. 
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We have considered that these modifications have to be small in order not to contradict 
the known tests of GR and, as a consequence, we have treated them as perturbations. 
Thus, we have applied a perturbative approach to compare the /(-R)-induced secular effects 
with the latest observationally determinations coming from various versions of the EPM 
planetary ephemerides. In particular, we have considered the ratios of the corrections to 
the standard secular precessions of the longitudes of perihelia estimated by E.V. Pitjeva for 
several pairs of planets in the Solar System. For all the models that we have considered 
(power law, Hooke-like force, logarithmic corrections, Yukawa-like force, gravito-magnetic 
effects) our results show that the perturbations deriving from the non-linearity of f{R) are 
not compatible with the currently available apsidal extra-precessions of the Solar System 
planets. Moreover, the hypothesis that the examined /(-R)-induced perturbations are zero, 
which cannot be tested by definition with our approach, is compatible with each perihelion 
extra-rate separately. 

This might suggest that, on the one hand, the /(i?)-induced secular effects cannot 
explain the observed extra-precessions and that, on the other hand, the effects of the 
non-linearity of the gravity Lagrangian are important on length scales much larger than the 
Solar System (e.g. on the cosmological scale) and their effects on local physics are probably 
negligible. It will be important to repeat such tests if and when other teams of astronomers 
will independently estimate their own corrections to the standard secular precessions of the 
perihelia. 
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